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ABSTRACT 



The ^^O/^^O ratio of the Solar System is 5.2 while that of the interstellar 
medium (ISM) and young stellar objects is ~4. This difference cannot be ex- 
plained by pollution of the Sun's natal molecular cloud by ^^0-rich supernova 
ejecta because (1) the necessary B-star progenitors hve longer than the duration 
of star formation in molecular clouds; (2) the delivery of ejecta gas is too ineffi- 
cient and the amount of dust in supernova ejecta is too small compared to the 
required pollution (2% of total mass or ~20% of oxygen); and (3) the predicted 
amounts of concomitant short-lived radionuclides (SLRs) conflicts with the abun- 
dances of ^^Al and "^^Ca in the early Solar System. Proposals for the introduction 
of ^^0-rich material must also be consistent with any explanation for the origin 
of the observed slope-one relationship between ^^O/^^O and ^^O/^^O in the high- 
temperature components of primitive meteorites. The difference in ^*0/^''0 ra- 
tios can be explained by enrichment of the ISM by the ^'''0-rich winds of asymp- 
totic giant branch (AGB) stars, the sequestration of comparatively ^^0-rich gas 
from star-forming regions into long-lived, low-mass stars, and a monotonic de- 
crease in the ^^O/^^O ratio of interstellar gas. At plausible rates of star formation 
and gas infall. Galactic chemical evolution does not follow a slope-one line in an 
three-isotope plot, but instead moves along a steeper trajectory towards an ^^O- 
rich state. Evolution of the ISM and star-forming gas by AGB winds also explains 
the difference in the carbon isotope ratios of the Solar System and ISM. 

Subject headings: astrochemistry — solar system: formation — ISM: evolution — 
Galaxy: evolution — stars: AGB and post- AGB 
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1. Was the protosolar cloud polluted by a supernova? 



Abundance ratios of the three stable isotopes of oxygen (^^0/^0/^0) in meteorites 
and planetary samples are used to infer the existence of reservoirs and processes in the 
early Solar System. Each of these isotopes was produced in a different mass-distribution 
of stars, and the isotopic composition of the Solar System is the cumulative result of ~9 
billion years (Gyr) of previous Galactic stellar nucleosynthesis. ^^O, a primary isotope, 
is principally formed during He-burning in massive stars, and added to the ISM by 
core-collapse supernovae (Type II SNe). ^^O and ^^O are secondary isotopes: The first is 
produced by reaction of light nuclei during the CNO cycle; the second by nitrogen reactions 
during He burning. Both of the heavier isotopes are ejected by Type II SNe. ^^O is 
also produced by the intermediate-mass progenitor s of asymptotic gi ant branch (AGB) 
stars. The yield scales with the abundance of ^^O (IMeyer et al. 20081 ). Thus, in a plot of 
^^O/^^O vs. ^^ O/^^O, Galactic che mical evolution (GCE) is expected to proceed along a 



slope-one line (ITimmes et al. 19951 ). Comparisons of the Sun with gas and other objects in 
the Milky Way should reveal such a trend, with younger objects being relatively ^^O poor. 

Observations to date do not fulfil this expectation. Measurements of primordial 
Solar System materials , including the solar wind returned by Genesis spacecraft 
(IMcKeegan et al. 20091 ). fall along a slope-one line that has a high ^^O /^^O (but not 



high ^^O) compared t o the much younger IS M and newly-formed stars (jWannier 1980l : 



Wilson fc Root 1994 



Wouterloot et al. 20081 ). One proposed explanation is that the 
Solar System formed from molecular gas polluted by -^ ^O-rich gas from m assive stars 



( lOlive fc Schramm 1982 



Henkel fc Mauersberger 1993 



Young et al. 20091 ). SLRs such as 



26, 



Al, Ca and Fe in the early S olar System are interpreted as evidence that contamination 



did occur (iGoswami et al. 20051 ). However, the oxygen isotopic composition of the Solar 



System falls along the same trend inferred for the initial compositions of AGB stars that 
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formed before the Sun and produced the presolar oxide grains in primitive meteorites 



( INittler 20091 ) . Furthermore, to explain the oxygen isotopic composition of the Solar 
System, (1) the delivery of SN ejecta must have occurred before cessation of star formation 
in the cloud, (2) it must have been efficient enough to produce the observed offset, (3) 
it must be quantitatively consistent with the abundance of SLRs, and (4) it must be 
compatible with the observed slope-one dispersion of primitive Solar System materials and 
the solar wind in a three-isotope plot (5^^0 vs. 5^^0), if this dispersion is primary (see 
below). We argue that the SN pollution hypothesis fails all four tests. 

First, stellar nucleosynthesis calculations show that only B-star progenitors (8- 
5 M(T,) have SN ej ecta with ^^0/^''0 ratios substantially higher than the current ISM 



( iMeyer et al. 20081 ) and are plausible sources of pollution. However, these stars have 



main-se quence lifetimes >15 Myr, longer than t he typical duration of star formation in 



clouds flWilhams et al. 2000 



Lada &: Lada 20031 ). and thus low- mass stars forming from the 



same cloud are unlikely to be polluted. The probability of an unaffiliated B star polluting a 
given cloud during its lifetime is approximately Rt f\ w here ~ 2 x 10^ Myr~^ is the Type 



II SN rate in the Milky Way (jPragicevich et al. 1999h fmo stlv from B star progenitors) 



r ~ 10 Myr is the mean cloud lifetime (ILada fc Lada 20031 ). and / ~ 10~® is the volume 
filling factor of a 10 pc cloud in the Gal actic disk. The probabil ity is ~0.2%, comparable to 
that of an encounter with an AGB star Jxastner fc Mvers I994I . 



Second, mixing between hot, tenuous SN ejecta and cooler, dens er molecular gas is 
inefficient. A 3-dimensional simulation produced mixing of 50 ppm (IBoss et al. 20081 ). 
only 2.5% of what is required. Refractory grains of AI2O3 can carr y ^^Al into a cloud 



(iGaidos et al. 20091 ). but can accommodate only 0.9% of all oxygen (ILodders 20031 ). SN are 
predicted t o be copious sourc es of dust, but no more than lO^^-lO""^ Mq has been detected 



per event (iMeikle et al. 20071 ). The scarcity of 0-rich presolar grains of plausible SN 
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origin (IHoppe &: Zinner 2000l ) also supports inefficient SN dust production and transfer to 
star-forming regions. 

Thir d, SLRs such as ^ ^Al ^^Ca, and ^^Fe, with half-hfes of 0.72, 0.1 Myr, and 



2.6 Myr (IRugel et al. 20091 ) . respectively, accompanied SN ejecta and Wolf-Rayet winds 



from massive stars, and S 



R abundances and any oxygen isotopic shift should be related 
( iGounelle &: Meibom 20071 ). Previous work has shown that introd uction of exogenous SLRs 
would likely have occurred during the giant molecular cloud phase (iGaidos et al. 20091 ). ^^Al 



was homogeneously 



flThrane et al. 2006 



distributed in the proto solar nebula with the ratio ^^A1/^''A1 ^ 5 x 10 
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Villeneuve et al. 20091 ). The allowed fraction of SN ejecta depends 



on the interval of free decay of ^^Al between injection and the formation of refractory 
calcium-aluminum-rich inclusions (CAIs, taken to represent time zero). The abundance of 
^^Ca constrains the interval of free decay to ~0.3 Myr, and Monte Carlo calculations show 
that the Solar System most likely formed from a clou d polluted to 0.3% by a generation of 



stars that formed 4.5 Myr earlier (IGaidos et al. 20091 ). Only progenitors more massive than 
50 Mq could pollute, but such stars are depleted in both ^''O and ^^O with respect to the 
Solar System. Limongi & Chieffi (2003) predict that the ejecta from a 35 Mq progenitor 
contains 5, 1.7 x 10~^, and 5 x 10"^ Mq of ^^O, ^''O, and ^^O, respectively. [See Woosley 
& Weaver (1995) for similar figures.] This ^^0-rich material leaves the ^^O/^^O ratio of the 
cloud essentially unchanged. Ejecta from higher-mass progenitors is expected to be even 
more depleted in ^^O and ^'^O. 

The high-temperature components (chondrules and CAIs) of un altered chondritic met e- 
orites plot along a slope-one hne in a three-isotope oxygen diagram ( Yurimoto et al. 20071 )11. 



^Only rare FUN (Fractionation and Unidentified Nuclear effects) and F (Fractionation) 
CAIs show evidence for extensive mass- depend ent isotope fract ionation due to melt evapo- 
ration and plot along lines with slope of 0.52 (IKrot et al. 20091 ) 
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The dispersion is usually explained by models in which (1) the initial oxygen isotopic 
compositions of both solids and gas were identical and equal to the solar wind value returned 
by Genesis; and (2) the slope-one line is a consequence of mixing with an ^®0-poor reservoir 
generated in the protoplanetary disk by "self-shielding" of CO isotopomers from ultraviolet 
dissociation, e.g. Lyons & Young (2005). However, these assumptions have been challenged 
by Krot et al. (2009), who concluded that (1) primordial (thermally unprocessed) solids 
in the Solar System were already ^^0-depleted relative to solar nebula gas; and (2) CO 
self-shielding had only a small effect on the isotopic chemistry of dust in the Solar System. 
In this scenario, the slope-one line was prod uced by mix ing of relict, isotopically distinct 
reservoirs, i.e. gas and dust of different ages (IDwek 20061 ). Any contamination by ^^0-rich 
ejecta would have to be contrived such that the final slope-one line was preserved. 



2. The Solar System in the context of Galactic chemical evolution 



The isotopic composition of the Solar System can be explained in terms of GCE that 
has deviated from the canonic al slope-one line at lea st since the Sun formed 4.6 Gyr ago. 



A flat age-metallicity re 
the solar neighborhood 



ation ( 



Holmberg et al. 20071 ) and the deficit of metal-poor stars in 



Caimmi 20081 ) indicate that the metallicity of star-forming gas has 



evolved little over the past ~10 Gyr, presumably because of the infall of metal-poor gas 
and the sequestration of metals in low- mass stars (iColavitti et al. 2008i) . This means that 
the recent evolution of the ^^O/^^O and ^^O/^^O ratios in the Milky Way depended not on 
^^O, but on the relative contributions by stars of different masses to each isotope, and the 
time-variation caused by the strong dependence of main sequence lifetime on stellar mass. 

The ISM is neither a closed system nor chemically and isotopically homogeneous. Gas 
and dust move through different phases (molecular clouds, neutral ISM, HII regions), and 
stars with different masses add or remove isotopes at different times depending on their 
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main sequence life. A GCE model should include: (1) infall of metal-poor, ^^0-rich gas; 
(2) massive progenitors of SN (and Wolf -Rayet stars) which begi n ejecting mass 1-3 Myr 



after the formation of their progenitors ([Limongi fc Chieffi 20031 ) and can pollute their 



natal molecular cloud (IGaidos et al. 20091 ): (3) stars with initial masses of 1.5-8 M© which 
move onto the red giant and asymptotic branches >70 Myr after their formation, long after 
dispersal of the molecular cloud, eject mass into low-density ISM phases, and preferentially 
enrich them in ^'^O with respect to star- forming region]^; and (4) stars with <1 Mq which 
remain on the main sequence for >12 Gyr and are a sink for all isotopes. The sequestration 
of relatively ^^0-rich gas from star-forming regions will leave the ISM enriched in ^^O. We 
propose that the cumulative effect is GCE with a stagnant ^^O/^^O ratio but a decreasing 
i8q^i7q j-atio, and that the difference between the Sun (and all stars of its age) and gas in 
the present Milky Way is a result of this gradual change. 

We calculated the evolution of the oxygen isotopic composition in the vicinity of the 
Sun's formation and present location using a two-box model of GCE. One box represents the 
low-density phases of the ISM (hereafter referred to as the ISM), which receives metal-poor 
infalling gas, gas ejected from star-forming regions, winds from AGB stars, and some SN 
ejecta. The second box represents star- forming regions (molecular clouds). It receives gas 
from the ISM, returning (most of) it after disruption of the clouds, and the difference is 
incorporated into stars. Three types of stars are considered: low-mass stars (< 1 Mq) which 
only sequester mass; intermediate-mass progenitors of AGB stars, and high-mass (> 8 Mq) 
SN progenitors. SNe eject mass into either molecular clouds or the ISM, depending on 
their main-sequence life relative to the cloud lifetime (10 Myr); AGB stars pollute only 
the ISM. The mass flux from the ISM into molecular clouds is calculated using a fixed 



One m anifestation of this is the appe arance of 0-rich presolar grains in primitive 



meteorites ( Nittler et al. 1997 



Nittler 20091 1. 
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ISM residence time of 280 Myr, a value calcula ted from a total gas s urface density of 
14 Mq pc~^ in th e present solar neighborhood (jNaab &: Ostriker 20061) . a star formation 



Lada &: Lada 20031 ). and a local star-formation 



efficiency of 10% flWilhams et al. 20001 : 
rate of 5 x 10~^ M^^ pc ~^ Myr~^. This last value is based on a SN r ate of 0.04 yr~^ 



( iDragicevich et al. 19991 ) and the abundance of Al in the 



This residence time is consistent with previou s estimates ( 



interstellar residence times of presolar grains (IHeck et al. 20091 ). The return flow of gas 



Milky Way ( 



Diehl et al. 2006h 



Tenorio-Tagle 20001 ) and the 



through HII regi ons and the disruption of molecular cloud s by SN assumes a cloud residence 



time of 10 Myr (iWilliams et al. 2000 



Lada Sz Lada 20031 ). We relat e star- formation r ate to 



total gas density by a Schmidt-Kennicut law with an index of 1.45 (IFuchs et al. 20091 ). We 
adopt a power-law stellar initial mass function with index of 2.35 for > 1 Mq and take the 
fraction of stellar mass in lower-mass stars to be 0.42. 

The present gas infall rate is poorly constrained by direct observations, but rates have 

been inferred from the distribution of stellar metallicities. We assume a constant rate based 

I I 
on a recent analysis of G- and K-star metallicities (iCasuso fc Beckman 20041) . Moreover, the 

adopted relationship between star-formation and gas s urface density, and the current total 

mass density in the solar neighborhood [~50 Mq pc~^ (jNaab fc Ostriker 20061 )] restrict the 

range of possible infall rates in our model. A rate of 3.5 x 10"'^ Mq pc~^ Myr~^ successfully 

reproduces the present surface densities of gas , stars, and stellar remnants [14, 35, and 3 



Mq pc ^, respectively flNaab fc Ostriker 20061 )]. The inf ailing gas is assumed to contain 



10% solar ^^O and no ^^O or ^^O. 

Nucleosynthetic yields are taken from Woosley fc Weaver (1995) for SN and from 
Boothroyd fc Sackmann (1999) for AGB stars. We do not consider stars more massive 
than 40 Mq, for which detailed yields have yet to be published. These stars are thought 
to contribute negligible ^^O or ^'^O. Dependence on metallicity is determined by linear 
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interpolations between solar and single sub-solar values. The ^^O yield from massive stars 
might have to be adjusted significantly downwards be cause of revised va lues for the proton 



capture reaction rates ^^0(p,a)^^N and ^^0(p,a)^^F (iMeyer et al. 20081 ). Removing the 
contribution from massive stars to ^^O altogether results in a model ^''O/^^O ratio in 
star-forming regions 4.6 Gyr ago that is within 15% of the solar value. With unadjusted 
^^O yields, the model predicts an ^^O/^^O ratio that is within 7% of the solar value. We are 
unable to simultaneously reproduce the surface densities of stars and gas, and the absolute 
oxygen abundance of star-forming clouds 4.6 Gyr ago (i.e. the solar abundance). One 
solution, besides tuning yields, is start the assembly of the disk at the Sun's location only 
10 Gyr ago, i.e. 3.5 Gyr after the formation of the Galactic center (INaab fc Ostriker 20061 ). 
Another is to allow a small fraction of the oxygen-enriched gas from star-forming regions to 
permanently escape the disk and be replaced with metal-poor gas. 

Figure [1] shows the predicted evolution of the oxygen isotope ratios in both the ISM 
(diamonds) and star-forming regions (squares) with respect to values in star-forming 
regions 4.6 Gyr ago (taken to represent the solar value). GCE is towards a ^^0-poor, 
but relatively ^^0-rich condition: it does not follow a slope-one line. In our model, the 
difference in ^^O/^^O between the ISM and star-forming regions is due only in small part 
to contamination of clouds by SN, because only progenitors of 18-40 Mq are sufficiently 
short-lived and these contribute little ^^O. Instead, it is mostly due to contamination of 
the ISM with ^^O from AGB stars over the ~280 Myr it takes for the two reservoirs to 
exchange. 

Episodic gas infall or star formation produces yet more dramatic departures from a 
slope-one trajectory. Cosmo logical models predic t stochastic time variation associated with 
infalling dark matter halos fiColavitti et al. 2008.) . and accompanied by elevated rates of 
star formation. If the Sun formed during a ~1 Gyr episode of increased star formation 
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(twice the background rate), it partakes of an enhanced ^^0-rich contribution from massive 
stars (Figure [2]) before winds from AGB stars return the ISM to its previous isotopic 
trajectory. In a sense, this is SN pollution, but on a Galactic scale. 



3. Discussion 

Our model of GCE in the solar neighborhood produces a super-unity slope in a 
three-oxygen isotope plot, and can explain the high ^^O/^^O rati o of the Sun with res pect 



to the ISM. The absence of this phenomenon in previous models (jPrantzos et al. 19961 ) may 
be attributable to the assumption that stars form from a gas of mean disk composition, 
rather than one to which AGB stars of the same generation have yet to contribute. Our 
model also predicts a difference between the con iposition of the ISM and star-forming 
regions: such a difference has not been observed (iYoung et al. 2009i) but measurements of 
oxygen isotopes are currently possible only in molecular (i.e. CO-containing) gas and the 
composition of the low-density phases of the ISM is unknown. The composition of young 
stellar objects also appears to be dispersed along a slope-one line in a three oxygen isotope 
plot. This effect could be produced by heterogeneous mixing of Galactic disk gas with 
infalling, metal-poor gas lacking ^''O and ^^O (Figures [Hand [2]). 

A purported Gal actocentric gradient in oxygen isotopes (increasing ^^O with radius) 



(IWilson fc Root 1994 ) is cited as support for standard G CE: However, recent data do 
not exhibit a significant trend (jPolehampton et al. 20051 ). and while the ^^O/^^O ratio at 



the Galactic ce nter may be higher tha n disk values, the former may have been previously 



overestimated (IWouterloot et al. 20081 ). Instead, Wouterlout et al. (2008) report an 
increasing ^^O/^'^O ratio with Galactocentric radius and a lower ^^0/^''0 in the Large 
Magellanic Cloud (LMC) compared to similarly metal-poor regions of the outer disk. They 



explain both of these observations as a result of the accumulation of ^''O with time, i.e 
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both the inner disk and LMC are older than the outer disk. A corollary of this conclusion is 
that at any given Galactocentric radius, younger objects (e.g., the ISM and newly- formed 
stars) will be more ^''0-rich than older ones (e.g., the Solar System). 

Our scenario predicts that most low mass stars of the Sun's age should have a similar 
oxygen isotopic composition. The SN contamination scenario, in contrast predicts that 
the Sun lies at the extreme of a distribution of ^^O/^^O ratios (corresponding to different 
levels of contamination), something the observations so far do not suggest. While it is 
presently not possible to measure t he oxygen isotopic comp osition of main-se quence stars. 



other isotopes such as those of Mg (lYong et al. 2004 ) or Si (ITsuji et al. 19941 ) might serve 
as proxies. 



The ratio of the Solar System (89±1) flWvckoff et al. 200d : Iciavton k Nittler 2004h 

is higher than the interstellar ratio in the solar neighborhood (76 ± 7) (IStahl et al. 20081 ). 
This difference can be explained by the high ^^C/^^C r atio in winds from red giant branch 



and AGB stars ( ITimmes et al. 1995 



in the ISM may also be important (IStahl et al. 20081 ). The present-day C and O isotopic 



Milam et al. 2005h . although chemical inhomogeneities 



composition of the ISM at the Galactocentric location of the Sun's birth is consistent with 
mixing of average AGB winds with the ISM 4.6 Gyr ago (represented by the Solar System) 
(Figure [3]). 

Our arguments do not preclude any pollution of the Sun's parent molecular cloud by 
massive stars, but demonstrate that the required amount and kind of pollution is neither 
plausible nor necessary to explain the ^^0/^''0 ratio of the Solar System. The Solar System's 
initial inventory of ^^Al from the ejecta of M > SOM© stars would have been accompanied 
by a shift of about -65%o in both heavy isotopes. If this ^'^O- and ^^0-depleted oxygen was 
delivered primarily as gas to the natal molecular cloud, it would have established an offset 
between the isotopic compositions of gas and dust in the cloud along the slope-one line. 
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This, as well as mixing with infalling metal-poor gas, c ould explain the s lope-one dispersion 



in primitive meteorite components and the solar wind (IKrot et al. 20091 ) . 



We thank Jonathan Williams and Ed Young for helpful discusssions, and Roberto 
Gallino for a detailed review. ANK and GRH acknowledge support from NASA grants 
NNX07AI81IG and NNX08AG58G, respectively. 
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Fig. 1. — Predicted evolution over 10 Gyr of the three oxygen isotope in the ISM (diamonds) 
and star-forming regions (squares) at 1 Gyr intervals, terminating in the upper right at the 
present day. Units are 1000 x the natural log of ratios with the value 4.6 Gyr ago when the 
Sun formed. The open ellipse is the composition of the ISM relative to the Solar System 



(IWouterloot et al. 20081 ). The dashed line is the standard slope-one GCE trajectory. The 
arrow is the mixing line produced by adding 20% gas with 10% solar ^^O and zero ^^O and 
^^O to present-day star- forming regions. 
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Fig. 2. — Same as Figure 1, except with a period of elevated (2x) star formation coinciding 
with the birth of the Sun. The event is modeled as a gaussian with a standard deviation of 
250 Myr. During the event the ■'^^0/-'^^0 ratio increases as more massive stars eject ^^O into 
the ISM; more slowly-evolving AGB stars add ^^0-rich gas later, returning the ISM and 
star-forming regions to their former paths. 
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Fig. 3.— ^^i-jQ 180/170 ratio in the Milky Way: The sohd hne is the Galacto- 

centric gradient in the ISM based on Milam et al. (2005) and Wouterlout et al. (2008). The 
circle is the Solar System value, and the compos ition of 1.5-9 Mq AGB star envelopes after 



second dredge-up ( iBoothroyd fc Sackmann 19991 ) are small squares. The large square is the 



IMF-integrated AGB ejecta composition assuming a Salpeter IMF and the ejecta mass func- 
tions of Boothroyd & Sackmann (1999). The dashed line is a hypothetical equal- metallicity 
mixing line between the mean AGB ejecta composition and interstellar gas 4.6 Gyr ago an d 



at a Galactocentric radius of 6.6 kpc, represented by the Solar System (jWielen et al. 19961 ). 
The mixing line passes close to the present composition of the ISM at 6.6 kpc (diamond), 
indicating that it can be produced by addition of AGB ejecta to the ISM at that location 
4.6 Gyr ago. 



